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We present a theory of the DC electron transport in insulators near Anderson-Mott tran- 
sitions under the influence of coexisting electron correlation and randomness. At sufficiently 
low temperatures, the DC electron transport in Anderson-Mott insulators is determined by 
the single-particle density of states (DOS) near the Fermi energy (Ep). Anderson insulators, 
caused by randomness, are characterized by a nonzero DOS at Ep. However, recently, the au- 
thors proposed that coexisting randomness and short-ranged interaction in insulators open a 
soft Hubbard gap in the DOS, and the DOS vanishes only at Ep. Based on the picture of the 
soft Hubbard gap, we derive a formula for the critical behavior for the temperature dependence 
of the DC resistivity. Comparisons of the present theory with experimental results of electro- 
static carrier doping into an organic conductor K-(BEDT-TTF)2Cu[N(CN)2]Br demonstrate 
the evidence for the present soft-Hubbard scaling. 

KEYWORDS: Mott transition, electron correlation, Anderson localization, randomness, disorder, Anderson- 
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Metal-insulator transition (MIT) has been one of the 
central issues in condensed matter physics. 1 ' When the 
electron correlation drives a MIT, e.g. as a Mott transi- 
tion, 2 ' it opens a gap in the single-particle DOS. On the 
other hand, Anderson transitions induced by randomness 
are characterized not by a vanishing carrier number but 
by a vanishing relaxation time. 3 ' Therefore, the Ander- 
son insulators exhibit no gap in the single-particle DOS. 

In real materials, however, electron correlation and 
randomness coexist inevitably. Let us consider how a 
single-particle gap behaves when randomness (e.g., spa- 
tially fluctuating potentials) is introduced into a Mott 
insulator. One might think that, when the disorder is suf- 
ficiently strong, localized impurity levels fill up and col- 
lapse the gap completely. This simple expectation, how- 
ever, is not correct; electron correlations open a soft gap 
in the Anderson insulator (see Fig. 1). 

A Coulomb gap proposed by Efros and Shklovskii 4 ' 5 ' 
(ES) is a typical example of such correlation-induced soft 
gaps. They showed that, assuming a nonzero DOS at the 
Ep, the ground state is unstable against a particle-hole 
excitation under the influence of the long-range Coulomb 
interaction. As a result, the DOS, A, vanishes toward Ep 
as 
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Fig. 1. (Color online) Schematic picture of single- particle DOS in 
Anderson-Mott insulator. 



repulsions and diagonal disorders within the Hartree- 
Fock approximation. 6,7 ' Even though only the short- 
range interaction is present in this model, we observed 
the formation of a soft Hubbard (SH) gap for d = 1,2, 
and 3. We found that the DOS is scaled in E near Ep as 



A(E)=aexp - {-7 log \E - E F \}° 



(2) 
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where E is energy and d is the spatial dimension. Since 
the ES mechanism is based on the excitonic effects, the 
Coulomb gap vanishes when the short-range part be- 
comes dominant, e.g., near MITs where the dielectric 
constant diverges. 

On the other hand, recently, the authors found an un- 
conventional soft gap with the help of numerical analyses 
on the Anderson-Hubbard model with coexisting on-site 
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where a, and 7 are positive constants. To clarify the 
origin of the SH gap, we assumed the existence of low- 
energy multiply-excited states with electronic structures 
that are globally relaxed from those of the ground state. 
Indeed, by considering the ground-state stability against 
these multiply excited states, we analytically reproduced 
the numerically observed DOS[eq. (2)] successfully. Such 
low-energy multiply-excited states are characteristic to 
random systems that have multivalley energy landscape 
or many (nearly-)degenerate metastable states. Indeed, 
the insulator with the SH gap is characterized by a 
nonzero spin-glass (Edwards-Anderson) order parame- 
ter, 8 ' 9 ' in contrast to the uncorrelated Anderson insu- 
lators; there exist many excited states nearly degener- 
ated with the ground state generically in spin/charge- 
glass phases. 

The long-range part of the Coulomb interaction is sup- 
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pressed to low energies near MITs. Therefore the effects 
of the short-range part practically determines the elec- 
tronic structures, and the SH mechanism becomes dom- 
inant in the experimentally accessible energy scale. The 
DC transport measurement is a useful and effective tool 
for investigating the DOS at low energies with a high res- 
olution in insulating phases. In the previous study, 7 ) we 
derived the temperature (T) dependence of the DC re- 
sistivity in the presence of the SH gap. Indeed, we found 
that experimental results for SrRui-^Ti^Oa are consis- 
tent with the SH scaling near the MIT. Recent develop- 
ment of electrostatic carrier doping using the field-effect 
transistor principle now allows us fine tuning of MITs 
in strongly correlated electronic systems. This motivates 
us to further investigate the DC electron transport in 
Andcrson-Mott insulators, in particular, the critical be- 
havior of the DC resistivity. 

In this letter, we clarify the electron transport near 
the Anderson-Mott transition in the presence of the 
SH gap. Particularly, we investigate the critical be- 
havior of the T dependence of the DC resistivity. We 
compare the present theory with a recent experimen- 
tal result of the electrostatic carrier doping into a 
thin-film single crystal of an organic superconductor k- 
(BEDT-TTF) 2 Cu[N(CN) 2 ]Br (, hereafter abbreviated as 
K-Br). 10 ) 

Figure 1 illustrates the single-particle DOS in an 
Andcrson-Mott insulator. For simplicity, we consider the 
electron-hole symmetric case here. High-energy extended 
states are separated from the low-energy localized states 
by a mobility edge located at the energy distance E g 
from Ep. The MIT is identified by a vanishing E g , which 
corresponds to a divergence of the localization length £ 
at E F . 

At high T i.e. k B T » E g , carriers are thermally ex- 
cited over the gap between the mobility edges, and they 
dominate the conduction. Thus the T dependence of the 
DC resistivity follows the Arrhenius law as 

PiT) = Poexp(^). (3) 

On the other hand, at low T compared to E g , the electron 
hopping between localized states near Ep becomes domi- 
nating because the density of the thermally excited carri- 
ers vanishes exponentially as T decreases for k B T < E g . 

Mott showed that, when there is a non- vanishing DOS 
at Ep, the conduction is dominated by hoppings between 
localized states within a narrow energy centered at Ep at 
sufficiently low T [variable-range hopping (VRH)] . In the 
VRH, the T dependence of the DC resistivity exhibits a 
universal behavior given by 



P(T) 



Po cxp 



k B T 



(4) 



Here Eq is an energy scale set by £ and A(Ep). Because 
the VRH explicitly depends on the DOS near Ep, the 
formation of the SH gap modifies the T dependence of 
the resistivity given by eq. (4) qualitatively. 

We start with a brief review of the phenomenological 
theory of the SH gap following refs. 6 and 7, and extend 



it for the later use. Now we consider a single-particle 
excitation from the ground state by adding an electron 
at a certain position/site. The SH theory claims that 
the excitation energy has a lower bound A (> 0) arising 
from the ground-state stability against multiply excited 
states. Note that, in uncorrelated Anderson insulators, 
A is always zero, leading to a non- vanishing DOS at Ep. 
Because A is determined by the energy scale of electron 
correlations between the added electron and surrounding 
electrons in the ground state, A is given by 



A exp (-6ft) . 



(5) 



Here, 6 oc and ft is the distance between the added 
electron and the electron nearest to it in the ground state. 
Within the single-particle picture, the wave function of 
an localized electron, <j>, is scaled in distance from the 
localization center r as 



r d/2 cx P (-r/a, 



(6) 



where £ is the localization length. Therefore, A is scaled 
in the mutual distance ft as 



(7) 



A oc | <j>{R) | 2 oc£- d exp(-2ft/£). 

By comparing cq. (5) and eq. (7), we obtain 

A oc C d - (8) 

On the other hand, the distribution of R is given by 

P(R) = a' exp (-b'R d ) , (9) 

where a' and b' are non-universal positive constants de- 
pending on the electron filling and £. By integrating the 
distribution of A with respect to energy, which is derived 
from eqs. (5) and (9), we obtain 



A(E) — a exp 



-7 log 



E-Ep 



(10) 



for \E — Ep\ < Aq. Here we take into account logarithmic 
corrections ignored in the previous studies, 6 ' 7 ) which is 
important in considering the electron transport near the 
MIT. Here a, and 7(= b' 1,>d /b) arc positive constants. 
Near the MIT, 7 is scaled with £ as 7 oc £, while a is 
non-critical. 

Following the discussion in ref. 7 and taking into ac- 
count logarithmic corrections omitted in the previous 
studies, we obtain the T dependence of the DC resis- 
tivity as 



p = po cxp < c 



exp [-c 1 |log(fc B T/Ao)| 1 / d ] 
k B T 



(11) 



where po, cq and c\ are positive constants. This SH scal- 
ing indicates that the DC resistivity diverges toward zero 
T slightly slower than the Arrhenius law; it reduces to 
the Arrhenius law when c\ = 0. Near the MIT, c and C\ 
are scaled with £ as 



Co 



Ao 



1 + 



0, (12) 



ci = d 1/d 7 _1 oc C 1 -> 0. (13) 
Note that, even near the MIT where the dielectric 
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constant diverges, the electron interactions remain long 
ranged at asymptotically long distances and low ener- 
gies. Thus we expect a crossover from the SH scaling to 
the ES scaling with decreasing T. Below the crossover 
temperature, the DC resistivity follows the ES scaling 4 ) 
as 

■ T X 1/2] 
l> = /'n<'X!> \^ 
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Fig. 2. (Color online) Arrhenius plot of experimental data for 
V e = 0, 40, 60, 80, 100, 120 V taken from Ref. 10. 



Recently, it was reported that a thin-film single crys- 
tal of an organic superconductor K-Br adsorbed to a Si 
substrate becomes a Mott insulator due to a negative 
pressure imposed by the incompressible Si substrate; 11 ' 
the substrate prevents the organic crystal from thermal 
shrinkage at low T. Subsequently, Kawasugi et al. re- 
ported successful doping of electrostatic carriers into the 
organic conductor by applying a gate-bias voltage. 10 ) 
Figure 2 shows the experimental data of the DC resistiv- 
ity at V g = 0, 40, 60, 80, 100, and 120 V taken from ref. 10. 
With increasing gate-bias voltage, the divergence of the 
DC resistivity at low T becomes weaker and weaker, in- 
dicating that the system approaches the MIT. This fine 
tuning provides us with a suitable system to test the 
present theory. 

Figure 3(a) shows an Arrhenius plot of the experimen- 
tal data for V g — V. The Arrhenius law fits well with 
the experimental data for T < 30 K at a first glance 
in the upper panel. However, on closer inspection in the 
lower panel, we find a clear deviation from the Arrhenius 
law; the data are convex upward for T < 40 K in the Ar- 
rhenius plot, indicating that the DC resistivity diverges 
toward low T in a rate slower than the Arrhenius law. 
Indeed, as demonstrated in Fig. 3(b), the data exhibits 
a perfect fit with the ES scaling rather than the Arrhe- 
nius law. This indicates a formation of localized states 
near Ep and the existence of a Coulomb gap under the 
influence of the long-ranged part of the Coulomb interac- 
tion. Thus the insulator at V g — V can be regarded as a 
Anderson- Mott insulator rather than a pure Mott insula- 
tor. The randomness in K-Br may originate from motions 
of the ethylene moiety in the BEDT-TTF molecule. 12 ~ 15) 
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Fig. 3. (Color online) Enlarged plot of Fig. 2 for V g = 0. The 
broken (black) line in (a) is a fit by the Arrhenius law, while 
that in (b) is a fit by the ES scaling. The lower panels in (a) and 
(b) show the experimental data divided by the fitted curves for 
the Arrhenius law and the ES scaling, respectively. 



Now we discuss the electron transport near the MIT, 
e.g., at V g > 0. Figure 4(a) shows an ES plot for 
V g = 0,40^60,80,100, and 120 V. As increasing V g i.e. 
approaching the MIT, the region where the data fit the 
ES scaling (ES region) is restricted to lower and lower T 
as shown by broken (black) lines for V g = 0, 40, 60, and 
80 V. Particularly, for V g > 100 V, the experimental data 
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Fig. 4. (Color online) (a) ES plot of experimental data (solid 
(red) curves) for V g = 0, 40, 60, 80, 100, 120 V. The broken 
(black) lines are fits by the ES scaling, and the dotted (blue) 
lines are fit by the SH scaling. The experimental data divided by 
the fitted curve for the SH scaling are shown for = 100, 120 
V in (b) and (c), respectively. 



become concave downward down to the lowest available 
T. This behavior is consistent with the expectation that 
the effects of the long-ranged part of the Coulomb inter- 
action are suppressed to lower and lower T near the MIT. 
Indeed, at a higher T above the ES region, as is shown by 
the dotted lines in Figs. 4(b) and (c), the SH scaling fits 
the experimental data much better for V g > 80 V. In the 
SH fit, it is numerically difficult to unambiguously deter- 
mine all the four parameters of the SH scaling. Therefore 
we assume A = 60, 50, and 40 K for V g = 80, 100, and 
120 V, respectively, so that Ao decreases with increasing 
V g . It should be noted that the fitting quality does not 
strongly depend on the choice of Aq. 



T T (\T\ 

Vg (v) 


Ao (K) 


Po 


log ( c ) 


Cl 


80 


60 


9.7 ±0.1 


4.8 ±0.1 


0.95 ±0.05 


100 


50 


10.1 ±0.02 


4.0 ±0.1 


0.95 ±0.02 


120 


40 


9±1 


4.6 ± 1.8 


1.7 ±0.8 



Table I. Fitting parameters for the SH scaling. 



The parameters of the SH scaling obtained by the fit- 
ting are summarized in Table I. The V g dependences of 
the fitting parameters (cq, ci) do not contradict the pre- 
dicted critical behavior [eqs. (12) and (13)], while the er- 
ror bars are too large to identify the critical exponents. 

In summary, we have investigated the DC electron 
transport in the insulators near the Anderson- Mott tran- 
sition. We have derived the formula for the critical be- 
havior of the T dependence of the DC resistivity based on 
the picture of the SH gap. The present theory was tested 
in detail against the recent experimental result on the 
electrostatic carrier doping into a thin-film single crystal 
of K-Br. Indeed, we have found that the experimental re- 
sult of the DC resistivity follows the SH scaling near the 
MIT, supporting the present theory. 

A direct observation of the Coulomb gap was reported 
by using electron tunneling spectroscopy in Si:B. 16 ) 
Electron-tunneling measurements allow us to further in- 
vestigate the critical behavior of the DOS, when they 
are combined with recently developed experimental tech- 
niques for the fine tuning of the MIT (e.g., the elec- 
trostatic carrier doping and X-ray irradiation in organic 
compounds 17-19 ) ) . 
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